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Abstract 

We consider the XXZ spin chain with diagonal boundary conditions in the 
framework of algebraic Bethe Ansatz. Using the explicit computation of the scalar 
products of Bethe states and a revisited version of the bulk inverse problem, we 
calculate the elementary building blocks for the correlation functions. In the limit 
of half-infinite chain, they are obtained as multiple integrals of usual functions, 
similar to the case of periodic boundary conditions. 
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1 Introduction 



Doped low-dimensional antiferromagnets have attracted a lot of studies especially since 
the discovery of high-Tc superconductivity. A particularly simple form of doping results 
from replacing some magnetic ions of the crystal by nonmagnetic one's. Open Heisen- 
berg quantum spin chains [1] are the archetype of one-dimensional models providing 
microscopic description of such systems. Indeed, the presence of non-magnetic impuri- 
ties into crystals having effective one-dimensional magnetic behavior has drastic effects 
on their low energy properties : the chain is cut into finite pieces with essentially free 
(open) boundaries leading to the breaking of translational invariance. As a consequence, 
physical quantities such as for example the magnetic susceptibility will get measurable 
corrections due to the presence of the boundary [2-11]. The same quantum spin chains 
have also acquired recently an important role in the study and the understanding of 
the interplay between quantum entanglement and quantum criticality [12-20]. There, 
the presence of boundaries also leads to noticeable effects, like in particular Friedel os- 
cillations [21-23] and the algebraic decrease of the boundary part of the entanglement 
entropy as a function of the distance to the boundary [19,20]. 

Correlation functions are central in the description of such effects and are in fact 
accessible in experiments. In particular the local magnetic susceptibility in the presence 
of boundary can be obtained using muon spin rotation/relaxation on the corresponding 
crystals, see for example [11]. More generally, correlation functions contain the neces- 
sary information to compare the microscopic models at hand to the reality, in particular 
through the measurements of dynamical structure factors accessible by neutron scatter- 
ing experiments [24-30]. While the computation of exact spectrum of Heisenberg chains 
has already a very long history, see e.g. [31-44], and references therein, computation of 
exact correlation functions of integrable lattice models such as Heisenberg spin chains, 
in particular out of their free fermion point where already considerable work was neces- 
sary [45 53], has been a major challenge for the last twenty years. Progress have been 
obtained using different routes and several results are now available for the correlation 
in the bulk, i.e., far from the boundaries [41,42,54-81], although still more progress 
is needed to obtain full answers. Advances have been obtained also in the presence 
of a boundary using in particular g-vertex operator methods [82, 83] and field theory 
approach [6-10,84-92]. 

The aim of the present paper is to develop a method to compute correlation func- 
tions of integrable open (finite and scmi-infinitc) spin chains in the framework of the 
(algebraic) Bcthc ansatz for boundary integrable models [93-105]. For this purpose, we 
will consider the example of the finite XX Z spin-1/2 Heisenberg chain with diagonal 
boundary conditions, including in particular non-zero boundary longitudinal magnetic 
fields, and its (semi-infinite) thermodynamic limit. Our results concern the general ele- 
mentary blocks of correlation functions at zero temperature, namely the average value 
of arbitrary products spin operators going from the boundary to an arbitrary site at dis- 
tance m from this boundary. Any correlation function can be written in terms of these 
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elementary blocks. Previous attempts towards this goal in the Bethe ansatz framework 
can be found in [7,106-108]. 

The strategy we will follow to solve this problem is closely related to the one used in 
the periodic case [59,60]. The central object in this approach is provided by the mon- 
odromy matrix of the open chain, which is a function of a complex spectral parameter 
A and of inhomogeneous parameters attached to each site of the chain. Following 
Sklyanin [95] , it is given as a quadratic expression in terms of the standard (bulk) mon- 
odromy matrix with the adjunction of the so called boundary K-matrix which encode 
the boundary conditions [95,96, 100]; in this paper only diagonal K-matrices will be 
considered. This boundary monodromy matrix satisfies a boundary Yang-Baxter alge- 
bra governed by two ii-matrices while the K-matrix itself satisfies its c-number version 
also called reflection equation [95] . These settings have been used by Sklyanin to extend 
the algebraic Bethe ansatz method to this open case. In particular, the Hamiltonian of 
the chain can then be reconstructed in terms of a weighted (with the K matrix) trace 
of this monodromy matrix. Hence as in the periodic case, one can consider a common 
set of eigenstates of the boundary transfer matrix and of the Hamiltonian. 

The first task towards the computation of the correlation functions is to identify the 
space of states of the open chain as generated by the action of the entries of the boundary 
monodromy matrix (depending on difi'erent spectral parameters Xj) on some reference 
state (here the state with all spins up or down); then eigenstates of the open chain 
are obtained from such actions thanks to the Bethe ansatz equations for the spectral 
parameters Xj [94,95]. Using this framework, we will show that it is possible to find 
determinant expressions for the scalar product between a boundary Bethe state and an 
arbitrary boundary state and consequently for the norm of the Bethe eigenstates. This 
is achieved along the lines used for the bulk case in [59] using the factorizing F-matrix 
basis [58]. 

The second problem is to obtain the action of the local spin operators on such states. 
In the bulk case, it was given by the resolution of the quantum inverse scattering prob- 
lem, namely, by the reconstruction of such local operators in site j in terms of a simple 
monodromy matrix elements (evaluated at A = multiplied from the right and from 
the left by products of the transfer matrices evaluated in the inhomogencity parameters 
for i = 1, J. In this bulk case, the Bethe eigenstates of the Heisenberg chain Hamil- 
tonian being also common eigenstates for the transfer matrix, it was straightforward to 
obtain the explicit action of the local spin operators on such Bethe states. 

The situation in the presence of boundaries turns out to be slightly more subtle : 
due to the breaking of translation invariance, boundary Bethe states are no longer eigen- 
states of the bulk transfer matrix, hence leading to a difficult combinatorial problem 
while using the expression of local operators described above. 
We solve this problem in three steps : 

(i) We first find a general (simple) relation relating boundary Bethe states to bulk one's. 

(ii) Then the reconstruction of local spin operators is obtained through a rewriting of 
the above quantum inverse scattering problem solution as a unique monomial in terms 
of the bulk monodromy matrix entries, avoiding in particular the presence of products 
of transfer matrices; it gives a new form for the general solution of the quantum inverse 
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scattering problem. 

(iii) Due to this new form of the solution of the quantum inverse scattering problem, 
the action of the local spin operators on any boundary Bethe state (expressed in terms 
of bulk one's) can then be given in a simple way and the result can be rewritten back 
in terms of sums of boundary Bethe states. 

Then, using scalar product and norm formulas for boundary states, we obtain any 
correlation functions as explicit sums of ratio of determinants of size half the length 
of the chain. In the thermodynamic limit (the limit of semi-infinite chain), these sums 
become multiple integrals with weights given in terms of the density of Bethe roots in the 
boundary ground state; this density function indeed describes the infinite size limit of the 
above ratios of determinants. For the so-called elementary blocks of correlation functions 
it gives proofs of the multiple integrals representations obtained previously [82,83] using 
the q-vertex operator method, here both in the massive and massless regimes of the 
chain. The problem of computing physical spin correlation functions will be addressed 
in a subsequent paper; this involves summing large number of the elementary blocks 
obtained here, using techniques similar to the one's developed for the bulk case [62,63]. 

This paper is organized as follows. In Section O we briefly describe the open XX Z 
chain with integrable diagonal boundary conditions and introduce the main notations. 
Section [3] contains some elementary algebraic properties of boundary operators and 
boundary states, and a description of the ground state in the thermodynamic limit. 
In Section [H the scalar products of Bethe eigenstates with arbitrary dual states are 
computed. In Section [5l we explain how, using a new version of the bulk inverse problem, 
one can derive the action of a product of elementary matrices on a boundary arbitrary 
state. Finally, in Section [6l elementary building blocks of the correlation functions are 
computed using the results of Section [5] and Section HI and we give their multiple integral 
representation in the thermodynamic limit. Some technical details are gathered in a set 
of appendices. 

2 The boundary XXZ chain: definitions and notations 

In this paper, we consider the XXZ Heisenberg spin-1/2 finite chain with diagonal 
boundary conditions. The Hamiltonian of a chain of M sites is given by 

M-l 

'^=Y1 {^mf^m+l+f^i^f^^+l+A((T^CT^+l-l)}+/l-af + /i+a|^. (2.1) 

m=l 

The local spin operators a^, am and at site m act as the corresponding Pauli 
matrices in the local quantum space "Km ~ C^, and as the identity operator elsewhere. 
The quantum space of states of the chain is !K = (Sim=i^m- In (|2.ip . A is the bulk 
anisotropy parameter, and h± denote the boundary fields. In what follows, they will be 
parametrized as A = coshry and h± = sinhry coth^±. 

To diagonalize the boundary Hamiltonian Ti., we use the modified version of the 
algebraic Bethe Ansatz proposed by Sklyanin in [95]. As in the case of periodic boundary 
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conditions, the eigenvectors are obtained as those of a family of commuting transfer 
matrices, which are constructed as fohows. 

Let i? : C ^ End(y F), F ~ C^, denote the i?-matrix of the XXZ model. 



= sinh(u + 77) i2(n), with R{u) 



/I 0\ 

h{u) c{u) 

c{u) h{u) 

Vo 1/ 



(2.2) 



in which 

b{u 



sinhn 



c{u) 



sinhr/ 



sinh(n + t?) ' ^ ' sinh(ti + r]) 
It is obtained as the trigonometric solution of the Yang-Baxter equatioE0, 
Ri2{u - v) Risiu - w) R23{v -w) = R23{v - w) Rniu - w) Ruiu - v). 



(2.3) 



(2.4) 



The i?-matrix satisfies the following initial, unitarity and crossing symmetry relations: 



i?(0) = V, 

Ruiu) R2ii-u) = 1, 
afR'^^{u-ri) a\ = -R2i{-u). 



(2.5) 
(2.6) 
(2.7) 



Here V is the permutation operator on V ®V , R21 = Ri2'Pi2, and *^ denotes the 
matrix transposition on the first space of the tensor product. 
Let also K{u;S,) be the boundary matrix 



K{u) = K{u;0 



sinh(n + ^) 

sinh(^ — u] 



corresponding to the diagonal solution of the boundary Yang-Baxter equation [96] 

Ri2{u - v) Ki{u) Ri2{u + v) K2{v) = K2{v) Ri2{u + v) Ki{u) Ri2{u - v). (2.9) 

A commuting family of transfer matrices T(A) E End3f is constructed from R and 
K as 



r(A) = tro{i^+(A)T(A)i^_(A)T(A)}. 



(2.10) 



Here the trace is taken over an auxiliary space Vq ~ C^, K±{X) = K{u it r]/2;^±) G 
EndVo, T{X) G End(Vo (8> 3i) is the bulk monodromy matrix. 



T(A) = RoAiiX - Cm) . . . i?02(A - 6) i?oi(A - 6), 



(2.11) 



^Here and in the following, indices label the spaces of the tensor product in which the corresponding 
operator acts non trivially. For example, in (|2.4p . which is an equation on Vi 55 V2 (81 V3 , Vi ^ C^, Rij 
denotes the i?- matrix (I2.2|l acting in Vii^Vj. 
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and T{X) is defined as 

f (A) = RioiX + ^i-v) ^2o(A + 6 -?/)•• • Rmo{X + - V)- (2-12) 

In these last expressions, Rom denotes the i?-matrix in End(Vb ® "Km), and .^i, ^2, • • • , Cm 
are arbitrary complex parameters (inhomogeneity parameters) attached to the different 
sites of the chain. Note that, due to (12. 6p and (12. Th . 

f(A)=7(A)ao^r*"(-A)ao^ (2.13) 
= 7(A)r-i(-A + ??), (2.14) 

with, in our normalization, 

M 

7(A) = (-1)^^ 7(A) = (-1)'^ n [ + ^^•) + - 2^)] • (^-l^) 

In the homogeneous limit (^m = ^?/2 for m = 1, . . . ,M), the Hamiltonian ()2.ip can be 
obtained as the following derivative of the transfer matrix (|2.10|) : 

2[sinhr/]i-2*^ d , , 

H = / / M r.^ / / M ^^(A)| + constant. (2.16) 

tr{i^+(r//2)} tr{K_(r//2)}dA ^1a=„/2 ^ ^ 

In the case of periodic boundary conditions, the space of states is constructed in 
terms of the operator entries A, B, C, D End IK of the bulk monodromy matrix 
(j2.1ip expressed as a 2 x 2 matrix acting on the auxiliary space: 



(A{X) B{X)\ 
^(^^-VC(A) DiX))- 



(2.17) 



These operators satisfy a quadratic algebra given by the following quadratic relation on 

Ri2{X - /i) Ti(A) T2{fi) = r2(/i) ri(A) i?i2(A - /x). (2.18) 

In this framework, eigenstates of the periodic Hamiltonian are constructed as the mul- 
tiple action of creation operators B{Xj) on the reference state |0) with all spins up, 
provided that the corresponding spectral parameters Xj satisfy the bulk Bethe equa- 
tions. 

In the case of the diagonal boundary conditions (|2.ip . a similar construction can be 
performed (see [95]) using the operators entries A-, C_, V- G End!K (respectively 
A+, B^, C-I-, P+) of one of the "double-row" monodromy matrices or Uj^ defined on 
End(Vb ® K) as 

U.{X) = r(A) K_(A) f (A) = (^^^Ix) vlU) ' (2.19) 
(A) = r*o (A) k!^ (A) f (A) = (^ll^J ^+ (^)^) . (2.20) 
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Note that the matrix W_ (respectively ^+), as weh as its operator entries A^, C_ 
and (respectively A-\-, B+, C+ and T>^) depend also on the parameters ^^i, . . . ,S,M 
and (respectively .^i, . . . ,£,m and ^_|_). It will be sometimes necessary in the paper 
to specify explicitly this dependency, denoting e.g. U-{X;S,-) instead of W_(A). The 
matrices and satisfy the boundary Yang-Baxter equations 

Ruiu - v) {U^)i{u) Ruiu + v-r]) {U-)2{v) 

= {U^)2{v) Ruiu + V-7]) {U^)i{u) Ruiu - v), (2.21) 
Ri2{-u + v) {U+)\' (n) Ri2{-u -v-v) {U+f^{v) 

= {U+Yi{v) Rui-u -v-rj) {U+)\' (n) i?i2(-n + v), (2.22) 

which leads to commutation relations for their operator entries. 

Note that the transfer matrices (j2.10p can be expressed either in terms of the matrix 
elements of U- , 

T(A) = tro{K+(A)Zi_(A)} 

= sinh(A + r?/2 + ^+)^_(A) - sinh(A + r?/2 - ^+) P_(A), (2.23) 

or in terms of the matrix elements of Uj^ , 

T(A)=tro{K_(A)Zi+(A)} 

= sinh(A - r]/2 + ^_) ^+(A) - sinh(A - r?/2 - ^_) P+(A), (2.24) 

and their common eigenstates can be constructed either in the form 

N N 



k=l k=l 



N N 



IV'-({A})) = n^-(^j)lO)' = (Oin^^Ai), (2.25) 

or in the form 

l^+({A})) = n^+(^J-)|0), (V'+dA})! = (OinC+(A,-), (2.26) 

k=l k=l 

provided the set of spectral parameters {A} satisfies the Bethe equations 

y,(A,;{A};e+,e-) = 2/,(-Ai;{A};e+,e-), J = l,...,iV, (2.27) 

with 

2/,(M,{A},C+,e-) - ,i,h(A,_/. + r?) sinh(A,+^-r?)' ^'"^^^ 
{A}; ^-) = -a(^) d{-fi) sinh(^ + ^+ - ??/2) sinh(/i + C_ - 7?/2) 

N 

X [sinh(;^ - Xk - v) sinh(/i + Xk - ??)] • (2.29) 
k=i 
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Here a(A) and d(A) stand respectively for the eigenvalue of the bulk operators ^(A) and 
D{\) on the reference state | 0), 



M M 



i(A) = J]sinh(A-ei + ??), d(A) = J]sinh(A-ei)- (2-30) 



i=l i=l 



The corresponding eigenvalue of the transfer matrix on an eigenstate ()2.25p or 

(12:26]) is 

' I sinh(2^)nL[KA.-/^)K-^-A,,)] 

+ a(-^)d(^)^^"^^'^ " ^^"^^^ " ^+ + ^/'^ ^^"^^^ " + ^/'^ 1. (2.31) 
sinh(2/.)nL[% + A,)6(M-A,)] J 

Let us finally introduce some convenient notations that we will use all along the 
paper: for any set of complex variables {xj}, we define 

Xjf^ = Xj X]^ and xjj^ ^ Xj ~\~ x^^. (2.32) 



3 Boundary states 

3.1 Algebraic elementary properties 

In this subsection, we collect some usefull elementary properties concerning boundary 
operators and boundary states. They mainly follow from the description of the boundary 
XXZ model in terms of the bulk one. Indeed, the "double-row" monodromy matrices 
l/{± of the boundary XXZ model being quadratic in terms of the bulk monodromy 
matrix T (see definitions (|2.19p - (|2.20p and formula (j2.13p ). the boundary operators are 
themselves quadratic in terms of the bulk operators. 

This quadratic nature influences non-trivially the dependence on the spectral param- 
eter of the boundary operators; in particular, a "Z2 invariance" arises in the spectral 
parameter dependence of the operators B± and C±. More precisely, the following propo- 
sition holds: 



Proposition 3.1 The boundary operators B± and C± satisfy the properties: 

^ smhSA+^g C_(-A) = -2^y|±^C_(A), 

smh(2A — rj) smh(2A — rj) 

B,(-A)=- -'°;;;^j;-"; B,(A). c4-A)=-2^ie^c,(A). 

smh(2A -|- rj) smh(2A -|- rj) 

Proof — Such properties are simple consequences of the boundary-bulk operator 
decompositions following from ()2.19p - (|2.20p . □ 
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Note that the proportionahty factors appearing in Proposition 13.11 are not intrinsic 
and could in principle be removed with an appropriate choice of the normalization of 
the iiT-matrix. 

This symmetry has important consequences since the operators B- or (respec- 
tively C- or C+) generate the quantum space of states of the boundary XX Z model 
by their multiple action on the reference state | ) (respectively, on the dual reference 
state (0 I). In particular, the previous proposition naturally suggests that the solutions 
of the Bethe equations (j2.27p are characterized by the same Z2 symmetry, which indeed 
can be shown from a direct study of the Bethe equations. 

Proposition 3.2 Let {Ai, . . . , Xjy} be a solution of the system of Bethe equations ()2.27p . 
then {cJiAi, . . . jCtatAat} is still a solution for aj = ±, j = 1, . . . , N . 

Proof — This follows directly from the form of the Bethe equations (12.27P . □ 

As in the bulk case, the operators entries of the boundary monodromy matrix can 
be related by some simple relations. This is the subject of the next lemma. 

Lemma 3.1 The following relations hold: 

<yo U±{X-, e±) = -T^U±{X- -e±) r,, (3.1) 

or explicitly: 

^±(A;e±) = -r:.P±(A;-e±)r,, C±(A;e±) = -r,S±(A;-e±)r,, (3.2) 

M 
fc=l 

Proof — These identities follow from the definitions (I2.19l) - (j2.20p and from the 
bulk identity T(A) = T(A) T^. □ 

The question now arises whether the state ( V'-({A}) | (|2.25p (respectively ( V'+({A}) | 
(I2.26P ) is actually related to the dual state of | ^_({A}) ) (respectively | ^/^+({A}) )), i.e. 
whether the operators B-{X) and C_(A) (respectively Bj^{X) and C+(A)) are conjugated 
to each other. Indeed, if the Hermitian conjugate V'^ of an operator V G End(Vo ® ^) 
is defined as 

V\X)=[V{X)Y'-'''\ (3.3) 

where *i --*a/ denotes the transposition on the quantum space 'K and * the complex 
conjugation on c-numbers, we have the following result. 

Proposition 3.3 In the vicinity of the homogeneous limit of the massless model (rj € 
iM, ^fc — 77/2 G M and (,± £ iW), U± (A) has the following Hermitian conjugate: 

ui{X) = -{U^{-X*)y\ (3.4) 

An analogous result holds in the vicinity of the homogeneous limit of the massive model 
(rj € M, — r//2 € iM and ^j- — iq±'K/2 € M with q± = Q^1 ), namely 

4(A) = (-1)^±{Z^±(A*)}*«. (3.5) 



where = ® 
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Proof — It follows from the conjugation properties for R and K: 

rI (A) = -<(-A*), KiX ± r?/2; e±)* = -K{-X* ± r?/2; 
in the massless case, and 

<(A) = <(A*), K{X±r,/2;^±r = {-l^^KiX* ±r,/2;U), 

in the massive case. □ 

In the next proposition, a set of formulae are derived to express the states of the 
boundary XX Z model in terms of those of the periodic bulk XX Z model. 

Proposition 3.4 Let Ai,...,Aiv be arbitrary complex numbers. Then the boundary 
states I ^£({A}) ) and {tps{{X}) \ , e = ±, can be expressed in terms of the bulk states as 

N 

(T1,...,0-]V=± j = l 

N 

{Mm\= E <\,...,a.)(^i'---'^iv;^-) (oin'^(^")' (3-7) 

(Ti,...,<TAr=± j=l 

where 

»,(^'. . . . .A^;?-) = n [ - 7(A,)<.(-AJ) «^ 



3=1 



■31 



l<r<s<iV 



^ i<r<.<iV smh(A^,j 



X 

^ J a 1 ./ /J 

l<r<s<Af 

N 



x...(AJ.^,-,/2)] n..^l!i|^. (3.") 



l<r<s<Af 

in which we have used the notations A^ = ajXj and X% = ajXj+akXk for j, k = 1, . . . , N . 
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Proof — Let us show (|3.7p and (|3.9|) for the state {ip-{{X}) \ by induction on N, 
the proofs for (■(/'+ ({A}) | and | ■0±({A}) ) being similar. 

For = 1, the expression fohows from the representation (IA.2P of the operator C_ 
and from the action of D on the dual reference state ( | . 

Let us now suppose that the decomposition (|3.7p - (|3.9p holds for any set of complex 
variables {Ai, . . . , Xn}- The action of C^{Xiy+i) on the state ( V'({Ai, . . . , Xn}) \ can be 
computed from (|A.2p using the expression of the bulk action of D(A/vr_|_i) [43]: 

N 

N N+l n sinh(Afcj + r]) 

( 1 n CiX,)D{XN+i) = diXk) ( I n ^(^^■)- (3.12) 

^=1 n sinhAfc,- J=l 

i=i 

In such a sum we can distinguish between the direct term k = + 1, and the in- 
direct terms k < A^ + 1. Let us show that each indirect term does not contribute. 
The indirect terms corresponding to a given A; < A^ + 1 are proportional to the states 
(0 I nf=ij^fc C{X]) C{Xn+i) C{-Xn+i) with coefhcients: 

s . , sinh(2Ajv+i - r?) ^ sinh(A^^ - + r?/2) 

7 Aat+i smhr/ — —7777 ^ > aN+id{akXk) . , ..^ — ^ 

smh(2A7v+i) smh(A^ - A^^^) 

sinh(A^ . + r]) . 

xff \ ^{ , fff- .(Ai,...,A^;^_). (3.13) 

-Li sinh(A^^.) ('^iv.'^iv)^ ^' ' Jv,s ; V J 

There we f actor ize the following expression: 

sinh(2AAr+i - rj) sinh(2Afc - ??) -pr sinh(Afc - AJ + 77) sinh(Afe + AJ - 77) 
sinh(2A7v+i) sinh(2Afc) sinh(Afc - AJ) sinh(Afe + AJ) 

X sinhry 7(A^+i) d{Xk) d{-Xk) ^f,^,,,...,,^)(Ai, • • • , A^v; e-), (3-14) 

which does not depend on the values of ak and (Ttv+i- Here, H,^^ • • • ? -^A^; 

is a part of H^~^ ^^-^lAi, . . . , Atv; <^_) which does not contain A^. The remaining sum 
in (|3.13p reads as 

^ sinh(A^ + - r?/2) sinh(A^^i - + r?/2) 

which is zero. 

Thus, only the direct action of ()3.12p contributes, and it generates (13. 7|) . □ 
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Remark 3.1 The above proposition implies that, for specific values of the spectral pa- 
rameters, the corresponding boundary and bulk states are proportional. For example, 
since d{^i) = a{^i — rf) = Q, we have. 



N 



|V^+(teJ)) = i7f+(teJ;e+) \lB{i^,m, (3.16) 

h=l 

N 

(V'-lteJ)! =i/f-(teJ;e-) (0| nctej, (3.17) 

h=l 

N 

|V'_(te, -r?})) =i/f-(te, -r?};e_) \{B{i,,^-r^m, (3.18) 

h=l 

N 

{^+m, - r/}) I = i?f+({e.. - (0 1 n - r?), (3.19) 

h=i 

in which {61, • • • , Ciiv) a subset 0/ {^1, . . . , ^a/}, and h'^^ [{\}]^,±) , for O = B,C, 
denotes the coefficient h'^^ i)(^1' • • • ^^N]i±)- 

Note that the previous proposition, together with the bulk decompositions of the 
boundary operators, allows us in principle to reformulate the quantum inverse problem 
for the boundary XXZ in terms of the periodic bulk one. Indeed, we will use this 
property, in Section [5l to compute the action of a product of local operators on a 
boundary state. 

Now, let us recall that the two expressions (j2.23p - (j2.24p of the boundary transfer 
matrix T coincide as well as the Bethe equations derived using the | V'- ) boundary 
states (I2.25P or the \ ones (I2.26p . In absence of degeneration, these observations 
naturally suggest that, for any solution of the boundary Bethe equations, the corre- 
sponding eigenstates | ■i/'-({A}) ) and | i/^+dA}) ) have to be proportional to each other. 
Indeed, this holds as shown explicitly in the next proposition. 

Proposition 3.5 Let {Ai, . . . , Aat} be a solution of the system of Bethe equations (j2.27p . 
Then the corresponding eigenstates generated by Bj^ and B- are proportional, as well as 
those generated by C+ and C_ .■ 

n 0) = n 7^1 _ OA n ^-i^m. (3.20) 

where 

r^ffx 1.^ „.^ _ A ^(Aj)sinh(Aj -x + r?/2) -p|- sinh(A^ + A, - 7?) 

G(|A4, X, y) - YY ^^^^^ ^^^^^^ ^ ^ _ ^^2^ ^^11^^ ^.^^^^^ ^ ^ . (3.22) 
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Proof — The above identities can be proved using the boundary-bulk decomposition 
of Proposition 13.41 by directly showing that the two ratios 

Hb+/B- - (-1) —sZ 77 r-' (3-23) 

do not depend on {cri, . . . , ctat} and coincide respectively with 

smh(2Aj+r/) w c ^ a TT smh{2Xj - r]) w c ^ 

^.^^ sinh(2A,-r,) H sinh(2A, + ,) GiM^i-^^^)' 

Let us consider for example the ratio Hq^/i3_ (|3.23p . which reads: 

^ _ sinh(2A,- + r?) rf(-AJ) sinh(AJ + g+ - r?/2) sinh(A^, - ??) 

n gi^i^(2A,- - r?) a(-Ap sinh(AJ - + r?/2) sinh(A-, + t?) ' 

The action of the transformation cJa ^ — ctq on such a ratio for a given a G {1, . . . , A^} 
gives: 

-V sinh(2A, + r?) d(Ag) sinh(A^ - e+ + r?/2)) ^ sinh(A^, + r?) 

'''^+/^- sinh(2Aa - v) a{K) sinh(A- + + r?/2)) sinh(A-, - r,) ' ^ ^ 

where -ffa,B+/i3_ is a part of H^^/^_ which does not contain A^. Proposition 13.21 now 
implies that {AJ, . . . , A^} is a solution of Bethe equations if {Ai, . . . , Xn} is a solution. 
Therefore, applying the Bethe equation (|2.27p for j = a, we have 

djXZ) sinh(A^ - g+ + v/2) ^ sinh(Ag, + r?) 
a(A-) sinh(A- + e_ - r//2) sinh(A-, - r,) 

_ d{-Xl) sinh(A^ + g+ - r?/2) -j^ sinh(Ag, - r?) 

a(-A-) sinh(A- - + r//2) ii sinh(A-, + r?) ' ^ ' ^ 

SO that (j3.25|) coincides with the expression of Hq^/q_. □ 



3.2 Description of the ground state 

The Bethe equations ()2.27p can be written in the logarithmic form as 

N 

2Mp{Xj) + 5(A,; C+, e-) + [^(^ifc) + ^(^ik)] = ^^rij, l<j<N (3.27) 

k=l 
k^3 
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where nj are integers (with Uj < rij^i), and where the momentum p, the scattering 
phase and the boundary contribution g are defined as 

. ^ sinh(A - g+ + r?/2) sinh(A - + t?/2) 

ff A;g+,g- = zln . ■ I c 7^- (3-30) 

smh(A + ^+ — r]/2) smh(A + 4_ — r]/2) 

In the homogeneous hmit, the corresponding eigenvalues of the Hamiltonian Ti in the 
spin M — N sector are 

AT 

^({A}) = sinh7y{coth^+ + coth^_} + 4 ^ { cosp(Aj) - A}. (3.31) 

In order to characterize the ground state of the half-infinite chain M ^ oo, one 
should distinguish the two domains — 1 < A < 1 (massless regime) and A > 1 (massive 
regime), for which we set: 

vr ~ vr 

«i = -^i' C = > 0, = with ~ 2 ^ - 2' for 1 < A < 1, 

aj = iXj, ( = -ri> 0, ^_ = -|_ + with f_ G M, for A > 1, 

in which 5 = 1 for < sinhC and (5 = otherwise. Thus, to a given set of roots 
{Xj} corresponds a set of variables {aj} given by the previous change of variables. 
Note that two sets of Bethe roots {Xj} and {ajXj}, where aj = ±, correspond to the 
same Bethe vector. Therefore, we consider only solutions {Xj} such that 5R(aj) > or 
3?(aj) = 0, 9(aj) < 0. 

The ground state of the half-infinite chain M — > oo has been studied in [102], [104]. 
It appears that the nature of the ground state rapidities depends on the value of the 
boundary field h- . 

In the case where ^_ < or ^_ > (^/2, the ground state of the Hamiltonian (j2.ip 
is given in both regimes by the maximum number N of roots Aj corresponding to real 
(positive) aj such that cosp(Aj) < A. In the thermodynamic limit M — > oo, these roots 
Aj form a dense distribution on an interval [0, A] of the real or imaginary axis. Their 
density 

p(A,)= lim [M(A,+i-A,)]-i (3.32) 

satifies the following integral equation: 
A 

p{X) + J [K{X - /i) + K{X + fi)] p{X) dX = (3.33) 
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Here, 

zvr 2tt smn(A + r/j smn(A — r]) 

and A = +00 in the massless regime, while A = — i7r/2 in the massive one. We may 
extend the definition of p as the solution of (j3.33p on the whole interval [—A, A]. It is 
then easy to see that p{—X) = /o(A). Therefore, p satisfies the equation 

A 

p{X) + j K{X- p) p{X) dA = (3.35) 

-A 

which means that the density of Bethe roots for the ground state of the open chain is 
twice the corresponding density in the periodic case. 

In the case < < C,, the ground state admits also a root A (corresponding to a 
complex a) which tends to r//2 — with exponentially small corrections in the large 
M limit. In that case, the real roots density is still given by (j3.35p . 

All the above results are valid in the homogeneous limit. However, for technical 
convenience, we also introduce a familly of inhomogeneous densities p(A,^), depending 
on an additional parameter ^, as solutions of the integral equation 

A 

p(A, i) + j K{X - p) p{p, dp = ^t(A, 0, (3.36) 

-A 

with 

sinh 77 

smh(A — 4j smh(A — q + f]) 
It is easy to see that p{—X, r] — ^) = p{X,^). Then the function 

1 

Ptot (A) = :ttt E K^' + P(^^ ^ - 



2M 

k=l 



satisfies the integral equations (j3.33p and (I3.35P in the inhomogeneous case, and tends 
to the ground state density p in the homogeneous limit. The inhomogeneous integral 
equation ()3.36p can be solved explicitely as 



P(A,6 



Csinhf(A-0 

°° /I ^2n\ 2 



in the massless regime. 



~ 11 — 2^ a~rr\ r\ — V' ? = ^ ' ™ massive regnne. 

^^=iVi + r / ^i(hA-6>9) 



Finally we would like to stress that all the functions in (|3.36p are holomorphic in a 
symmetric strip of width r] around the interval [—A, A]. Therefore this equation still 
holds at the extra root A when it exists. 
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4 Scalar products of boundary states 



In order to compute correlation functions following the method proposed in [59], [60], it 
is necessary to have an explicit expression of the scalar products between a Bethe state 
and a general state. In the bulk case, such scalar products have been represented in the 
form of a determinant of usual functions in [109], [59]. The method proposed in [59] has 
been used in [107] to obtain a similar representation for the open XXX chain. In this 
section, we give the explicit expressions of the scalar products of boundary states in the 
XX Z case, and explain briefly how to derive them. 

4.1 Partition functions 

It is useful, as a starting point for the computation of scalar products, to consider the 
following functions: 

(M, {6}; e±) = ( I B±{Xm) . ■ ■ S±(Ai) I ), (4.1) 
2&({Aa},{6};e±) = (0|C±(Am)...C±(Ai)|0), (4.2) 

where | ) is the reference state with all spin down and ( | is its dual. Note that these 
functions correspond to the partition functions of the six- vertex model with domain wall 
boundary conditions and one reflecting end. 

Proposition 4.1 The above partition functions are related to each others in the fol- 
lowing way: 

({Aa},fe};^±) = (-1)*'2&({AJ,{6};-C±), (4.3) 
z!;{{Xa},m;C+) = (-1)^^47 ({-A«},fe};^+), (4.4) 
4tiM,m;^+) = (-i)^^^ff ({-AJ,fe};^+). (4.5) 

Proof — As Fa; I ) = I ) and ( | Fa, = ( | , the relations ()4.3p are direct conse- 
quences of Lemma [3. 11 The other two identities can be proved using the boundary-bulk 
decomposition of Proposition 13. 4[ □ 



Proposition 4.2 [110] The partition function can he represented as the determi- 
nant 

M 



2m({A4, {&};?-) = n [7(A/3)a(A^)a(-A^)] 

/3=1 

M M 

n n [ sinh(A/3 - Cfc) sinh(A^ + Cfc)] 

f3=lk=l 

Yl [sinhA/j^sinhA/3^] fl [ sinh sinh(^^^, — rj) 

/3<7 r<s 



detAAC-(A„,efc;e-), (4.6) 
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where 

j^- ^ sinhr? sinh(2Aa - r/) sinh(^^ + jk - ri/2) 

"'^ smh(Aa -ik + ri) sinh(AQ + Cfe - ??) sinh(Aa - ik) sinh(AQ + ^k)' 

Proof — In [110], both a set of recursion relations and the corresponding solutions 
were obtained. 

We propose in Appendix[B]an alternate derivation of this representation. It is based, 
just as in [59], on direct calculations in the basis induced by the twist F introduced 
in [58]. □ 



4.2 Scalar products 

Let us define, for two sets of complex variables {Ai,...,AAr} and {^i, . . . , ^utv}, the 
following different scalar products: 

ST' ({A}; {/^}) = ( ({A}) I v.. m) ), (4.8) 

for El, 62 € { + , -}. 

Theorem 4.1 Let {Ai, . . . , A^r} he a solution of the system of Bethe equations (2.21) 
and {/Ui, . . . , /iAf} he generic complex numhers. Then, the scalar products between the 
state \ip+{{n})) and the eigenstate {ijj-{{X}) \, and hetween the state {ip+{{fi})\ and 
the eigenstate \ , are respectively given as 



5^'^({A}; M) = n b(Aa) d{Xa) d(-A.)] ^j^^'^ ^||^|; ||;|| , (4.9) 

4'"(M;{A}) = nb(^'^)«(^'^)«(-^'^)] tt'^^vliAiiM) ' (^-'^^ 

where the matrices 7 and V are defined as 



T«/3({A},M) = ^t(/.^,{A}) (4.11) 

V (f\lSnl\- sinh(2A„) sinh(2^^ - r?) 

Vapli^f, ifJ^i) - gi^j^(2A„ - 7?) sinh(/i^ - K) sinh(^;j + A„) ' ^ ^ 

in which T{iip, {A}) denotes the eigenvalue (]2.3ip of the transfer matrix T(/x) on a Bethe 
eigenstate parametrized by {Ai, . . . , Aat}. 



Proof — Let us for example consider the scalar product (j4.9p . This formula can 
be proved following the same procedure as in [59] for the bulk case. As the reference 
state is invariant under the action of the operator F, we can rewrite this scalar product 
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in the F-basis and use the exphcit representations C_ and given in Lemma 17.11 for 
the operators C_ and jB+ in this basis: 

N N 

5V+({A};M) = (0| nC-(A„)n^+(^'')|0)- (4-13) 



a=l 6=1 



The idea is then to insert, in front of each operator a sum over the complete set 
of spin states . . . , i^a ); where . . . , ) is the state with ni spins down in the 
sites ii,. . . ,im and M — m spins up in the other sites. We are thus led to consider 
intermediate functions of the form 



s+1, ■■■ ,^N) 



N 



im+1, ■ ■ ■ ,iN)- (4.14) 

a=l 6=1 

which satisfy the following simple recursion relation: 

, /Xi, . . . ,fj,m,im+l, ■■■ ,iN) 

= X] {j,im+l,---,'iN\B+{flm)\im+l,---,'i'N) 

X G('""^)({Aa:},//i,... ,/Um-i,j, >^Af)> (4-15) 

Note that the last of this function is precisely the scalar product we want to compute, 
G(^)({Afc}, ....fiN) = ( V'-({A}) I V'+IM) ), (4.16) 
whereas the first one, 

N 

GW({Afc},ii,...,^Jv) = (OinC-(Aa)|ii,...,iiv>, 

a=l 

is closely related to the partition function computed in the previous section: 

GW({Afc},ii,...,i^)= n ] X{HK-ii)h{-K-ii)\^h-\i,^-ii)\ 

l^ii,...,iM a=l /3=1 J 

X zJ-({Ai, . . . , Xn}, fe, • • (4.17) 

Solving the recursion (j4.15p we obtain, when particularizing the result to the case 
m = N, 

N 

N N sinh^r? n [7(Aa)d(A,)d(-A,)] 



a=l 



''=1 n [sinhAab sinhAab sinh/Uba sinh/lj 

a>b 

^ sinh(2Ab - rj) sinh(2^b + ^) ^ ^ u ^/l18^ 
X n det/7.,({A}, {/.}), (4.18) 
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with 



rr (fxx f \\ 7(/^/3){?/a(/^/3; W) " Vai-f^p; {A})} 

sinh(AQ - Hf^) smh(Aa + fip) 



This ends the proof of (14. 9p . 

As for (|4.10p . it can be proved following the same procedure, provided one writes 
the operators C+ and in the F-basis as in Lemma l7. 11 □ 

All the remaining scalar products can be expressed in terms of those given in The- 
orem [4?T1 Indeed, we have the following corollary. 

Corollary 4.1 Let {Ai, . . . , A^v} be a solution of the system of Bethe equations (2.21) 
and {/ii, . . . , //Tv} be generic complex numbers. Then, 

5V+(M; {A}) = 5-+({-A}; {-^}), (4.20) 
5^''({A}; M) = {-A}). (4.21) 

Finally, the proportionality between =b Bethe eigenstates, given in Proposition \3.5l al- 
lows us to complete the list of scalar products where one of the boundary states is an 
eigenstate. 

Proof — The idea is to go to the F-basis, and then to insert, between the product 
of the operators C and the identity as a sum over convenient intermediate spin states, 
and finally to use the results of Proposition 14.11 

For example, applying this procedure to the left hand side of ()4.20p in the F-basis, 
one obtains the relation 

N N N N 

( 1 n c- {^^b) n ^+(^«) 1 ) = ( 1 n ^-(-Aa; ?+) n ; i « (4-22) 

h=l a=l a=l 6=1 

which holds for any arbitrary sets of complex numbers {A^} and {//fe}. The Bethe 
equations and the scalar product formula (|4.9p being invariant under the simultaneous 
exchanges — > ^- and ^_ we obtain, under the additionnal assumption that 

{Aa} is a solution of the system of Bethe equations and thanks to Proposition 13.21 that 
the scalar product on the right hand side of (j4.22p is given by (j4.9p evaluated at {— Aa} 
and {— /ifc}. □ 



4.3 The Gaudin formula and the orthogonahty of Bethe states 

The scalar product formulae derived in the previous section can be used, as in the bulk 
case, to compute the norm of boundary Bethe states and prove the orthogonality of 
eigenstates corresponding to different solutions of the Bethe equations. 
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Corollary 4.2 Let {Ai,...,AAr} and {fii, . . . , fi^} be two different solutions of the 
system of the Bethe equations, that is, 



{cJiAi, . . . ,0-ArAAr} 7^ {/xi, . . . ,/XAr} for each aj = ±, j = I, . . . , N. (4.23) 

Then, for ei,e2 € {+, — }> the scalar product S^'^^ {{X}; {fj,}) vanishes, i.e. the corre- 
sponding Bethe states {ip£-^{{fi}) \ and \ ipe2{{X})) are orthogonal. 

Proof — Let us show the orthogonahty of the Bethe states corresponding to two 
different solutions {Ai, . . . , Xn} and {/ii, . . . , ^in} of (|2.27p . In such a case, the scalar 
products S^''^^{{X}; {fj,}) are proportional to each others for the different values of ei, £2, 
and the orthogonality follows from the fact that the determinant in (j4.9[) is equal to zero. 
Indeed, there exists a non-trivial vector ^({A}, {fi}) such that: 

N 

Y,Hjk{{X},{fi})vki{X},{f,}) =0, for any j = l,...,iV. (4.24) 

k=l 

Such a vector f ({A}, {/i}) can be constructed in the following way: 

N N 

n sinh(Aj - Hk) n sinh(Aj + fik) 

Vj{{X}, {n}) = — — = . (4.25) 

n smhXjk n smhAjfc 

To check that the equations (j4.24p are satisfied with the vector ()4.25p , we use the explicit 
expression for the matrix elements Hjk{{X}, {/J-}) (j4.19p and apply the Bethe equations 
for the solution {/ii, . . . , /iAr}. □ 



Corollary 4.3 Let {Ai, . . . , Xn} be a solution of the system of Bethe equations. We 
have 

N 

n b\K)diXa)di-Xa)ya{-Xa;{X})] 

S^'^m-, {A}) = sinh^,?^=^ 

n [ sinh Xab sinh Xab] 

^ ^ sinh(2A,-,) sinh(2A, + ,) (4^6) 
j}^ sinh2(2Aj) ^ j/cvi 

Here <!>' is the Gaudin matrix: 

^' (S\\\ ^ 1 j/fc(-Afc;{A}) 

^ikKW) = log , (4.27) 

dXj yk{Xk;{X}) 

with yj{x;{X}) defined as in Ii2.29\) . The norm of the corresponding Bethe eigenstate 
follows then from Proposition \3.5\ and Proposition [57 
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4.4 Partial scalar products in the thermodynamic limit 

For the computation of correlation functions, it is usefull to have an expression for partial 
renormalized scalar products. For some partition U a_ of |1 ; A^], we consider the 
sets of variables {Ai, . . . , Xn} and {/ii, . . . , /iAr}, with {A} solution of the system of Bethe 
equations (|2.27p . such that 

{A} = {Aa},ea_ U {Xb},^^^ , {^} = Maea^ U ' (4-28) 

in which {^i^}bea+ is a subset of {^i, . . . ,^m}- Then, 

'g^'^({A};M) ^ TT 7(Ab) sinh(2Afe) sinh(2Afe - r?) sinh(26, + r?) ^(g^,; {A}) 
5^'^({A}; {A}) 7(e^J sinh(2^ij sinh(2^i, - r?) sinh(2Afe + r/) y(Ab; {A}) 

^ -TT sinh(Aba) sinh(Aba) -i-r sinh(AfcQ) sinh(Afea) detAf 
a,b<L+ ^^^^'^^ibia) sinh(4ij sinh(^i, - Xa) sinh(^i, + Xa) det^v AT ' 

a>b fee«+ 

with y given by (|2.29p and 

Mab = 2M6ab{p{Xa) + ^giXa, C+ , " ^ E [^(^« " Afc) + ^'(Aa + Afc)] 

fc=l 

27r 1 

+ — i^(2Aa)| + 27r[K(A„ - Xb) - K{Xa + A^)] , (4.29) 
Mab=l^'''' if &Ga_, ^^^^^^ 

It remains to caracterize the ration of the two determinants of A4 and Af, wich reduces 
to the determinant of a matrix S of size |a+|: 

It is actually possible, as in the bulk case, to compute explicitely Sab for the ground 
state in the thermodynamic limit. 

Indeed, it is easy to see that, if Xj corresponds to a real root aj, 

Op real 

This follows from the fact that, if A^ corresponds to a real root of the ground state, the 
matrix element Afjp can be expressed as 

M,p = 2M6,p{p{Xj) + + 27r[i^(A,- - Xp) - K{X, + A^)], (4.33) 
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from the symmetry property p{X, fi) = p{—X, r] — p) of the inhomogeneous density, and 
from the inhomogeneous integral equation (j3.36p . 

Therefore, if the ground state contains only real roots (i.e. in the case ^_ < or 

L > c/2), 

'^"^ aJ^oo 2M^) • (^-^^^ 

Let us now consider the ground state in the case < ^_ < C/2, and let Ai = A 
corresponding to the complex root (i.e. Ai — > rj/2 — ^_); then, there exists 7 > 

M^oo 

such that 

Mip = g'{X){6ip[l + 0{Me-^^')] +0(e-^*^)}, (4.35) 

with [s''(A)]~"'^ ~ isinh(A — ^_ + rj/2) in the large M limit. Therefore, using (I4.32p . the 
inhomogeneous integral equation (j3.36p at the point A and the estimation of {Af~^)ai 
following from (I4.35p . we obtain that 

N 

Sab = {J^-')al Mlb + (-^"')a/3 M/3b, (4.36) 

ITT sinh (A - C_ + 77/2) [p(A, J - p(A, rj - J] if a = 1, 



if a / 1. 



2Mp{Xa) 

5 Action of local operators on a boundary state 

In order to compute correlation function, one should now determine the action of the 
corresponding local operators on a boundary state. As in the bulk case [59], the idea is 
to solve the inverse scattering problem, i.e. to express local operators in terms of the 
generators of the Yang-Baxter algebra. 

A natural idea would be to try to express these local operators directly in terms 
of the generators A+, C+, (or B-, C_, V^) of the boundary Yang-Baxter 
algebra. However, although it is quite easy to reconstruct in such a way a local spin 
operator at the first site of the chain [106], it seems much more difficult, due to the lack 
of translation invariance, to obtain effective formulas on the other sites of the chain. In 
pratice, the reconstruction proposed in [106], which involves the adjoint action of the 
bulk translation operators {A + -D)(Cfc), is unadapted to compute correlation functions 
of the boundary model since eigenstates of the Hamiltonian are no more eigenstates of 
these translation operators. 

Quite surprisingly, it is actually possible to use directly a version of the bulk inverse 
problem to compute the action of local operators on a boundary state. In this section, 
we will show how to reformulate the bulk inverse problem so as to circumvent the use of 
the translational invariance of the chain. Then, using Proposition 13.41 we will act with 
the corresponding products of bulk operators on a boundary state, obtaining a sum over 
some bulk states that eventually reduces to a sum over boundary states. 
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5.1 The bulk inverse problem revisited 

Let us define a familly of algebra homomorphisms ^ ^ ~* ^ as 

Xi : X h->- . . . RiiXRii . . . Ri-ii. (5.1) 

Then, for a local operator Xi at site i (i.e. which acts non trivially only on !Kj), Xi{Xi) 
can be expressed in terms of the bulk monodromy matrix entries as 

X^{Xi) = [a{ii) d{ii - 7?)] tro [To(eO ^o] {A + D){ii - r?), (5.2) 
= H^^) d(^^ - r,)] {A + Dm tro [a^' - X^] . (5.3) 

To compute the bulk correlation functions, the authors of [59], [60] used the reconstruc- 
tion (j5.2|] FI. together with the fact that 

i-l i-l 

X, = \{{A + D){i^)x^{X,) \[[{A + D){^^)]-\ (5.4) 

This was convenient there because the product of bulk transfer matrices merely produces 
a numerical factor when applied on a bulk Bethe state. As it is no longer the case when 
applied on a boundary Bethe state, the strategy is to simplify this product instead. 
We can make the following observation: 

Lemma 5.1 The products of the bulk monodromy matrix elements T^, f,i {^i) Tg gi {^i — rj) 
and Te'eiii - v) Tg' ei^i) vanish if e = e. 

Proof — It follows directly from the fact that Xii^i ^) Xii^i^ ) = ^/x./i Xii^i ^ ) a-^d 
from the expressions (|5.2|) . (j5.3|) of Xi (-^f > Xi {^i^ ) terms of the bulk monodromy 
matrix, in which ^ denotes the elementary matrix at site i with elements (^E^ ^ )ab~ 
^a,fi ■ □ 

Remark 5.1 Other interesting identities may be proved in the same way. For example, 
from the fact that Xi{EP) Xi{Ef^) = Xi{E}^) Xi{El^) , one obtains that C{Ci)B{Ci-v) = 
-A{ii)D{ii-r^). 

This result can be generalized to a product of 2m operator entries of the bulk 
monodromy matrix: 

Theorem 5.1 For any set of inhomogeneity parameters {^ii, ■ ■ ■ jCin}; the following 
product of bulk operators 

Te,„ iiir.) ■ ■ ■ Te^, (^n) (^^i " ^) • • • T,^^ e-,^ (Ci™ " ^) (5-5) 

vanishes if, for some A; G {ii, . . . , Ek = £k- 

^Note that we express here the resuh in a slightly different form, using that [(A + D)(^i)] ^ = 
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Proof — It can be proved by recursion on n, the case n = 1 following from Lemma [5. II 
Let us suppose that the result holds for n — 1, and consider the product (j5.5p with 
Ei^ = Ei^ . Using the commutation relation given by the quadratic algebra ()2.18p , one can 
move the exterior operators (at position n) T^. Tg {^i^ — rj) through those 

at position n — 1 (evaluated respectively at and —f])- Considering all possible 
cases for (^j ), T^. (fj — ?7)andT,. ifi .),T/.. (£j r?),itis 

easy to see that the resulting product of the 2(n — 1) inner operators should vanish due 
to the recursion hypothesis. □ 

This leads to the following corollary concerning the reconstruction of a product of 
local operators acting on successive sites of the chain: 

Corollary 5.1 A product of elementary matrices acting on the first m sites of the 
chain can he expressed as the following product of entries of the hulk monodromy matrix: 

m 
i=l 

X (6) • • • T,,^ [U) Te^ iU-v)--- Te, (6 " v) (5.6) 
with Ei = e[ + I (mod 2). 

Proof — This is a direct consequence of the solution (15. 2p , ()5.4p of the inverse problem, 
of the fact that [{A + D){C,)]~^ = [a{Ci) d{^i-r])]~^ {A + D){(i-r]), and of the previous 
theorem. □ 



5.2 Action on a bulk state 

Let us now establish the action of a product of elementary matrices of the form (j5.6p 
on an arbitrary bulk state | {Aj}i<j<Af ) = Ilj^i 0)- We refer for example to [59] 

for the explicit expression, in our notations, of the action on such a state of a single 
operator entry of the monodromy matrisH. Let us just recall that, like in (j3.12p . the 
action of A(^) or -D(/i) produces two kinds of terms: a direct term, which leaves the 
state untouched, and indirect terms, resulting in new states with one Xj replaced by fi. 
With this terminology, the action of operators of the form (|5.6p can be computed using 
the following lemma: 

Lemma 5.2 The action on a hulk state \ {Xj}i<j<M ) of a string of operators 



with El = e'i + 1 (mod 2), has the following properties. 

^There the action on the left was considered, but the coefficients are the same when one considers 
an action on the right instead. 
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• The only non-zero contributions of the tails operators (2) come from 

(i) the indirect action of all A{^i — rj) operators; 
(a) the direct action of all D[^i — r]) operators. 

• In what concerns the head operators (1), 

(Hi) if e'l = 1, the action of the operator T^^'^eiiCi) (i-^- ^Hi) or B{^i)) does not 
result in any substitution of a parameter — t]; 

(iv) if e'l = 2, the action of the operator T^i^^^{£^i) (i.e. D{£^i) orC{S,i)) substitutes 
— V with S^i; moreover, if there were others parameters — rj, j ^ I, in the 
initial state, they are still present in the resulting state. 

Proof — (i) and (Hi) follow from the fact that a{^i — 77) = 0. 

In order to prove (ii), let us consider the action of some operator -D(^i, — rj): its 
indirect contribution produces a state of the type -Bfe; —ri)\ {lJ'j}i<j<N-i ), for a certain 
set of parameters {fJ-j}. However, Theorem 15. II guarantees that the operator product 

is zero, hence (ii). 

Let us now prove (iv) by induction. If e[_^ = 2, then T^^ ^.^(Cii) acts on a state of 
the form -B(^ji — r/)| {/^j} ) for a certain set of parameters {fJ-j}, and is either equal to 

• -D(^j^), which acts only indirectly, and which can only replace — rj with ^j^; 
indeed, any other replacement would produce a state of the form — 
77)1 {fij} ) (where {p-j} is a subset of {fJ-j}), which is zero according to Lemma [5?T1 

• (^j), which gives, using Remark |5.H 

CteJ ^te, - V)\ {/^,} ) = -A{^^,)D{^,, - r?)| ); (5.8) 

D{^i^ — rj) acts only directly since A(^i-^^)B{(^i-^ — rj) = from Lemma 15. H and 
A((^i-^) cannot replace any other — t] since a(^j — rj) = 0. 

If (iv) is proved for all operators until position / — 1, then, from (i), (ii), (Hi) and the 
induction hypothesis, T^,/ ^. ((^jj acts on a state of the form B{(^i^ —rj) \ {/Xj} ) for a certain 

set of parameters {^j}, and the same reasoning as for I = 1 applies. □ 

In order to express the action of a product of elementary operators ITjli-^/' ^ 
a bulk state, let us consider the following set of indices: 

f3+ = {j : I <j < m, Ej = 1}, card(/3+) = s', 
l3^ = {j : 1 < J < m, e'j = 2}, card(/3_) = s. 
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Since our final goal is to compute correlation functions, one considers here only products 
such that the total number of indices in the sets /3+ and /?_ is s + s' = m, as otherwise 
the corresponding ground state average value is zero. We can thus introduce a set of 
indices ip G {1, . . . , m} such that 

/3- = {ip}p(z{i^,„^s}, with ik < ih ^or < k < h < s, (5.9) 
/?+ = {ip}pg{s+i,...,m}, with ik > ih ioT s < k < h <m. (5.10) 

From Corollarv 15.11 and Lemma 15.21 we obtain the following result: 

Proposition 5.1 The action of a product of elementary operators on an arbitrary bulk 
state can be written as 

m , N N+m 

l[Epl[B{Xk)\0) = J2^Pmi^i,---^N+m) n B{\k)\0), (5.11) 

j=l k=l f3,n k=l 

in which we have defined Xn+j ■= S,m+i-j for j € {1, . . . ,m}. In 115. the sums are 
over all the sets of m indices (3rn = {bi, • • • , bm}, where the bp are defined by 



(5.12) 



\bp G {l,...,N}\{bi,...,bp^i} forO<p<s, 
|6p e {1, . . . , N + m + 1 - ip} \ {bi, . . . ,6p_i} for s < p < m, 

and the coefficient J-p^ is 

N 

m r X n sinh(Afeb. +T/) jv . ^ . 

k=l ^ 
k^bj 

m 

. , . s n sinh(A,,^ -S.k + r]) 

smhiXb b +V) ^ ■ u^x ^ ^ 

i<i<j<m ^ ^ p=i n smh(A6j, - ^fc) 

m 

n sinh(Cfc - Afep + r/) 

X n • (5-13) 

p=s+i Yl smh(^fe - Xbp) 

Let us point out that, if the parameters Ai, . . . , Aat are solutions of the bulk Bethe 
equations, such a result agrees with what can be obtained with the method used in r^n'RI 



''taking into account that we consider here an action to the right, whereas in [60] we considered an 
action to the left. 
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5.3 Action on a boundary state 

We use the decomposition ()3.6p . (|3.1ip of boundary states into bulk ones in order to 
compute the action of a string of elementary operators on an arbitrary boundary state 
constructed from operators. It is remarkable that we are eventually able to express 
explicitely the result as a linear combination of such boundary states. Indeed, using the 
same notations as in Proposition 15.11 we have, 

Proposition 5.2 The action of a product of elementary operators on a boundary state 
can be written as: 

m , N N+m 

J{^PJ[B+{\um = Y.T;j{\]) n ^+(Afe)|0), (5.14) 

j=l k=l I3m k=l 

k^0m 

with Prn defined as in (|5.12p and the coefficient given as 

^ ^a{\l) ({\^\) ^ sinhAf. 



X 



,1 CLiij) H?+(if I) sinh(Af. +7?) 

n-i-r f -TT sinh(AJ + eXj - rj) „ sinh(^j + eXi) 



, . sinh(A![ + eXi) sinh(f,- + eXi — n) 

n sinh(AJ^ + ??) 

^ TT f TT sinh(gj - AH iea+ 

.11 1 11 sinh(Cj - Af - r/) n sinh(A^.) 

m m 

^ n smh{Xl-^k + v) m n sinh(a. - A^^ + r/) 

nk=ip + l k=ip + l fr 
^ li 7^ • (5-15) 

P=i sinh(A^^ - 6) P=«+i sinh(a-A^^) 

k^N+m+l-hp 

Here, the sum is performed over all aj G {+, — } for j € a+, we have defined Af := crjAj 
for i S Pm, with cij = 1 if i > N , and 

a+ = [3m^{l,---,N], = {1, . . . , iV} \ a+, 

^_ = {N + m + l- i}je/3„n{7V+i,...,iV+m}, 7+ = {1, • • • > "i} \ 7_. 

The function H^^ {{X}) is the coefficient (j3.1ip appearing in the boundary-bulk decom- 
position. 

6 Elementary building blocks of correlation functions 
6.1 Finite chain 

It is now a matter of straightforward calculations to derive the expression of elementary 
building blocks of correlation functions at zero temperature. They are given as the 
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ground state average value of products of elementary operators of the form 



^fr^ ^- (^.({A})i^.({A})) ' ^'-'^ 

in which Ai, . . . , Xn denote the ground state rapiditie^. 

Using Proposition 15.21 and the partial scalar product expression of Section 14.41 we 
obtain for the finite chain: 



Proposition 6.1 The boundary elementary building blocks can be written as 

n sinh^H n sinh(Cp<7 - ??) 



m TV N N+m N+m H+^{{X}) 



m which 



(-1)^' n n n smh(A- + - ^) „ . ^ , 

iJ+ fr\i)-V^ '=1^=1 ^ -TT sinh(gfc + - r?/2) 

{fe.}^^ Yl sinhiXl,^ + r?) sinh(A^^,^, - r?) J-^i sinh(A- + - r?/2) 



J l<i<j<m 

X 



if f "P ^ Hi >. 

n n - n ^^"^(^^p - ^ + ^) 
p=i fc=i fc=ip+i ^ 

X J] I J] sinh(A^^ - 4) n sinh(A^^ - Cfc - r?) [ detJ7. (6.3) 

p=s+l ^ k=l k=ip+l ' 

In this expression, the sum is performed over all € {+, — } for bj < N, and abj = 1 
for bj > N, and the mxm matrix Q is given in terms of the matrix S of Section \4^\ as 

^ik = -5N+m+i-bi,k, for bi > N, (6.4) 
^ik = Sb„k, forbi<N. (6.5) 

6.2 Half-infinite chain 

Let us now consider the thermodynamic limit M ^ oo of this quantity. 

In the case where all the roots aj describing the ground state are real, i.e. if ^- > C/2 
or ^_ < (see Section l3.2p . the sums over the indices bj from 1 to become, as in 
the bulk case, integrals over the density of the ground state. More precisely, we have to 
perform the replacement 

^ A A 

E / dA.^(A.) E ^.■/(\^) = / dA,/(A,)/,(A,). 



M 



"^-^ -A 



^Note that, due to Proposition 13.51 we could have also chosen to compute this average value by 
means of states - instead of states +. 
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Moreover, the sums over bj > N can be written as contour integrals thanks to the 
identity 

2z7rRes/9(A,0| =-2. (6.6) 

In the region < 2^_ < ^, one should also take into account the existence of the complex 
root. The term of the sum which corresponds to A can also be written as a contour 
integral since 



Res 



-sinh(A + ^_ -ri/2) 
Therefore, one obtains: 



I . = !• (6.7) 



(-1)"'-^ 



n sinh(^ij) n sinh{^ij - rj) 

j<i i<j 



s 

X 



s „ m 

n dA, J n dA,- Hrr^{{Xj} ; {6}) det [$ (A,-, a)] , (6.8) 



i=i ^ j=s+i 



with 



and 



^ (Ai,a) = I [p(A,-,^fe) - p{Xj,ri- a)] , (6.9) 



n n sinh(Aj + a - ??) ™ ■ ^rt , t /o^ 

o- ^ r ^ T . i ^ ^ i=l ^=1 TT smh(a: + - r?/2) 

l<i<j<m 
X 



W I [J sinh(Ap - ik) n - Cfe + ^) [ 

^ n 1 n - n ^^^^^^^ - - [ • (6.10) 

p=s+\ ^ fe=l fe=ip+l ^ 

The contours C and C depend on the boundary field. They are defined as 

f[-A,A]ur(A) ifo<e-<C/2, 

1 [ — A , A] otherwise , 

c = cur({4}). (6.12) 

where r(A) (respectively Td^;;.})) surrounds A (respectively ^i,...,^m) with index 1, 
all other poles being outside. 
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Remark 6.1 One can easily verify, as a consistency check for the above formula, that 
the reduction property 

{ El''''^ . . . E'-''-e];^_^, ) + { El^'''^ . . . E'-^''-e'J_^, ) = ( El^'''^ . . . E'-^^'- ) 
from m + 1 sites to m sites is satisfied. 



Let us finally rewrite explicitly this result in the two different regimes (massive and 
massless) of the XX Z model, using the fact that, in both regimes, the determinant of 
the matrix $ can be calculated explicitly (the corresponding expressions are given in 
Appendix ICj) . 

In the massless case, one gets directly: 

m _ 

n cosh(^ea) n [sinhd^z) sinh(fe"fcO] . 
^ fi ' ' n sinh(C^,) n sinh(e,, + iC) / V C 

j<i i<j C 

m ,. mm •l/\lA,■A^ 

n/oAj \ -pr -pr Smh(Aa + gfc + <) 

,=,+1 ^ K ^ M H si^l^ f (^'^ - ^k) sinh |(Aa + Cfc) 



X 

c 



sinh(|A/fc) sinh(|A,fc) ^ sinh(pfc) sinh(.efc + 4 + ^- 



sinh(AH - iQ sinh(AH + iQ sinh(Afc + 4 + 

J] I Yl sinh(Ap - 6) n sinh(Ap - 6 - [ 

p=l ^ k=l k=ip+l ' 

m f ip — ^ m s 

J] I J] sinh(Ap - Cfc) n ^"'^(^ - ^fc + <) [ , (6.13) 

)=s+l fc=l A;=ip+1 ^ 

in which C = CU T{{^k}), with 

C=^ , fo,«_<0or«_>C/2, 
lMur(-i(C/2 + e-)) forO<C_<C/2. 



X 
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In the homogeneous limit = —i(/2, these elementary blocks have the following form: 

(n4''i)=(-ir---%=^Binh".j_(^)"''"^7ni^. / n ^ 

sinh (|Afci) sinh [fX^i] ^ sinh (|Afc) 



X 



fe<i 

n 



sinh(Afei - <) sinh(Aw + iQ j^J^ sinh(Ajk + 4 + C-) 
' sinh"*+^''-i (Ap + i|) sinh'"-^'' (Ap - 4) 



cosh^- (f Ap) 

sinh-+v-^ (Ap + 4) sinh-v(Ap + zf) 
Ul cosh- (f A,) • ^'-''^ 

To obtain the explicit expression of the elementary blocks in the massive regime, one 
performs the change of variables aj = iXj, = i^k- Hence, using the corresponding 
representations for the determinants of the matrix one obtains: 

m _ 

n [O3{Pa)04iPa)] U. [OliPkl) OliPkl)] 
^£j,e, a=l k<l 



r * 

X 



/ TT [ TT ('^\ TT TT siii(«a + h + K) 



i\aik) 9i{aik) fV Oi{ak)02{ak) sin (^^ + 4 + ^^-) 



sin(aM - iQ sm{aki + iC) ^^^^ sin (a^ + i| + 



X n j n ^^^i^p - (^k) n ^^^^^^ -Pk- o ^ 

p=l fe=l k=ip+l ^ 

m ^ ip — ^ m >. 

X J] I []sin(ap-/3fc) J] sin(ap - A + [, (6.16) 

p=s+l ^ k=l k=ip+l ^ 

in which ^^(A) = 6i{X, q), with q = . The integration contours are C = CLI r{{Pk}), 
with 

^ ^ f [-7r/2, 7r/2] for |_ < or |_ > C/2, ^ 

- \[-7r/2,7r/2]ur(-z(C/2 + f_)) for < < CA 

In the homogenous limit Pj = —i(/2, the elementary building blocks for the correlation 
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functions are given as 



YIe'/''' ) = 2™(™+l) Sinh™^- q-^^^^ JJ [(1 + g^")^™ (1 - g2n)m(3m+l) 
j=l n=l 

r^^daj\ f /doA 0i{aki) Oijaki) 

x(-l) yiU^^jy .11 l2^J li sin(aH - ^C) sin(aH 



11 sin (afc + 4+ i^.) 11 
- sin-+v-i ^ .0 ^.^^-.^ (^^ ^ 

p=s+l 



-4 



Let us finally remark that all these computations can also be performed in the case of 
an external magnetic field along the /S^ direction. In that case, the integration contours 
and the density function will depend, like in the bulk case, on this external magnetic 
field. 
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7 Appendices 

A Boundary creation and annihilation operators 

Using the quadratic relations (j2.19p - (j2.20p . one can express the boundary operators 
A±, B±,C±,'D± in terms of the bulk operators. Note that it may sometimes be more 
convenient to rewrite the creation and annihilation boundary operators B± and C± in 
the form 

B.{X) = -7(A) [B{-X) A{X) sinh(A + C_ - ^/2) 

+ S(A)A(-A) sinh(A-e-+r//2)], (A.l) 
7(A) sinh(A + C_ - r//2) 

+ Z)(A) C(-A) sinh(A - + r//2)] , (A.2) 



C-(A) = 7(A) ''''^I'^^^^^A)'^^ [^^"-^^ ^^^^ ''""^^^ + " "^/^^ 
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and 



^+(A) = 7(A) "'"'^^l^^'l!^^ [Bi-X) D{X) sinh(A - e+ + r//2) 



+ B(A) D(-A) sinh(A + e+ - W2)] , (A.3) 

C+(A) = -7(A) + [yl(-A) C(A) sinh(A - e+ + r?/2) 

sinn(2Aj 

+ A(A) C(-A) sinh(A + e+ - r//2)] . (A.4) 

It is then convenient, for the computation of partition functions and scalar products, 
to express the boundary operators in the F and -F-basis. The concept of factorizing F- 
matrices was defined in [58], following the concept of twists introduced by Drinfel'd in 
the theory of Quantum Groups [111], and we refer to [58] for the explicit construction of 
the F and i^-matrices in the periodic XX Z spin-1/2 chain and for the representations 
of the bulk operators in the F and F-basis. Using the F and F-basis expression of the 
bulk operator, one obtains the following result: 

Lemma 7.1 Let X± denote the expressions of the boundary operators X± in the F-basis, 
and X± their expressions in the F-basis. Then, 

^+(A) = -f;n(-A,e.|e+)ar5 (^(-^-^^^(^-^^■) « ) , (A.5) 
i=i V u V [j] 

and 

- ^ /I A 

B_(A) = - g „(A, « _ _ . (A.7) 

where 

sinh ri sinh(2A - r?) sinh(x + £,-r]/2) 

n(A,4 j;) = 7(A) a(A) a(-A) —7-. z , . , : . (A. 9) 

^ ' ^ /w w V ; smh(A - ^ r/) smh(A ^ - ry) 



B Partition function 

In this appendix, we propose a proof of Proposition 14. 2[ Similarly as in [59], this 
derivation is based on direct calculations in the basis induced by the twist F introduced 
in [58]. Indeed, in this particular basis (called F-basis), the explicit expressions of 
the bulk operators A, B, C, D simplify drastically. Since the boundary creation and 
annihilation operators are quadratic in terms of the bulk operators, they have themselves 
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much simpler expressions in this -F-basis (see Appendix 1X1 for details). As moreover the 
states ( I and | ) are respectively invariant under the left-action of F and the right- 
action of -F~^, the partition function can be directly written in the -F-basis as 

2m ({Aa}, m-, e-) = ( I C_ (Ai) . . . C_(Am) I ), (B.IO) 

where C_(A) = FC_(A) F'^ is the expression of C_ (A) in the basis induced by F. Using 
now the expression (|A.6p of the operator C_(A), one obtains a new recursion formula for 
the partition function, which corresponds to a development of the determinant in (j4.6p . 
Indeed, acting with C_(Aj\/) on the state | 0), one has 

M 

47({Aa},Ui};e-) = ^n(AM,e*|C-)(0|C_(Ai)...C_(AM-i)IO, (B.ll) 

i=l 

where | i ) is the vector with all spins down except in site i and where the expression 

of u{Xm , is given by formula (|A.9|) . Since (crf)'^ = 0, the action of the other 

operators C_(Aa), 1 < a < M — 1, on the vector | i ) is diagonal on the space i, so that 

C— 

we obtain the following recursion formula for ZjJ : 

M 

i=l 

X Zlj_^{{Xa}a^M,{Cj}j^uC-)- (B.12) 

The coefficient of the recursion is 

M-l 

CM{XM,^i,{^j};^-)=u{^M,^i\^-) ll[b{Xk-^i)b{-Xk-^i)] llb-\C,i), {B.13) 

k=l jy^i 

which, as a meromorphic function of Xm, can be rewritten as 

cm{Xm, Ci, {Cj}; = 7(Am) a(AAf) a(-AA/) sinh(AM - sinh(AA/ + Ci) 
^ ^j-r^ sinh(A/3 - smhjXp + yr sinh(AM - ^j) sinh(AAf + ^j) 
sinh Am/3 sinh Am/3 fj{ sinh sinh(^jj - r]) 

M-l 

<f;,(AM,e.;e-)- 2^5/3 A/J-(A/3,e.;e-) . (B.14) 



with 



_ 1 / sinh(2A/3 + t]) sinh(2A/3 - rj) 

sinh 2A/3 sinh(2A/3 - r/) \ sinhAM/3 smhXMis 

n [ sinh AAffc sinh AAffc] m . c \ ■ u(\ 

k=i -TT smh(A/3 - 0) smh(A^ + ^i, 

X igri 11; 



k=l 
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This actually corresponds to the development with respect to the last line of the deter- 
minant 

M 

({Aa},{C,};e-) = n [7(A/3)a(A^)a(-A;3)] 

3=1 
M M 

n n [ sinh(A;3 - 6) sinh(A^ + 6)] 

/3=lfc=l 



n [ sinh A/3^ sinh Xp^] \[ [ sinh sinh(^s^ - r/)] M 



detAf^-, (B.15) 



r<s 



where J\f^- is the matrix obtained from Af^- by substracting to the last line Lm the 
linear combination of the other lines Yl!^=i Q^^p- Thus, as A/"*^" and J^^- have the 
same determinant, this concludes the proof. 

C Determinant of the densities 

In this Appendix, we give the explicit expression of the determinant of the matrix $ 
involving the density function of the ground state. 
In the massless case it is: 



rn 

det [$(A,-,a)] = [-J"\[[sinh(^^Xa) cosh 



a=l 



c 



n [sinh(|4z) sinh(|^fe;) sinh(^A;fc) sinh (|Aifc)] 

X '^' mm , (C.16) 

n n [ sinh f(A„-efe) sinh f(A„ + Cifc)] 

a=l k=l 

Let us now compute the determinant of the densities in the massive case, where the 
density of Bethe roots can be written in terms of theta functions: 



n>l 



(C.17) 



with q = = e ^. We therefore have to compute the following determinant: 



det [«>(A,-,efe)] 



' n>l ^ 



2m 



det 

m 



02{aj-f3k) , 02{(yj+Pk) 



+ 



(C.18) 



with aj = iPj, Pk = i^k- 

Let us consider detm [$(Aj,^jfc)] as a certain function / of the variable ai. It is an 
elliptic function of order 4m with periods tt and 2iC,. An irreducible set of poles is 



(C.19) 
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whereas 



±02, . . . , ittm, ±"2 +<,•••, iOm + < (C.20) 

are zeros of /. / being an odd function, A = is also a zero and, since / (A) = 
— / (A + iC), so is A = Up to congruence, there remain two other zeros which differ 
by iC) say xq and xq + i^. Since the sum of the zeros is congruent to the sum of the 
poles, Xq is either congruent to or to 7r/2. In fact, the only choice compatible with 
the periods of / is j;o = 7r/2. This means that, up to a constant independent of ai, f 
can be f actor ized as 

m 

n [9i{au)ei{au)] 

ei{a,) 02(ai) . (C.21) 

n [ei{ai-Pi)ei{ai+pi)] 

i=l 

This argument can be easily extended to all ckj, 1 < j < m, thanks to the antisym- 
metry in these variables. We can also apply a similar procedure to the variables fik, the 
difference being that we now deal with an even function and that the extra zeros are 
— iC/2 and — '7r/2. Finally we obtain 



^ lit 

det [^(A,,^)] = (--)" n 

i=l 



• (C.22) 

n [9i{ai-P,)e^{ai + pj)] 
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